Comments on "Length scale dependence of DNA mechanical properties 
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In a recent paper by Noy and Golestanian [l[ (NG) 
the elastic properties of DNA were studied by molecu- 
lar dynamics (MD) simulations. Two important conclu- 
sions were made: (i) the computed bending and twisting 
rigidities of DNA agree well with experimental data, and 
(ii) for lengths of a few helical turns, DNA dynamics 
exhibits long-range correlations in qualitative difference 
from the worm-like rod (WLR) model. Earlier similar 
studies showed that (i) the current MD forcefields sys- 
tematically overestimate the DNA rigidity, and (ii) MD 
trajectories of DNA involve only short-range correlations; 
no deviations from the WLR model are detectable if MD 
data are analyzed properly. Here it is argued that the 
data analysis in NG was not correct and that the earlier 
conclusions are valid. 

The first note concerns ergodicity of MD sampling and 
statistical errors. This is a standard issue for such stud- 
ies, but it was omitted in NG. In the most intensive earlier 
simulations the DNA rigidity was evaluated from 164 ns 
trajectories of one helical turn 0, 0] ■ This duration ex- 
ceeded the measured relaxation times by about 3 orders 
of magnitude and provided the accuracy below 5% for 
the computed persistence lengths. NG analyzed 100 ns 
dynamics of 3 and 5 DNA turns and it was implied that 
by considering many internal stretches of long DNA one 
improves the sampling. However, this intuitive assertion 
is valid only for stretching. For bending and twisting 
the true gain in sampling is close to zero or even nega- 
tive because the relaxation of a double helical fragment 
is slower when it is a stretch of longer DNA [1, [Bj] . In ad- 
dition, the self-diffusion coefficient of SPCE water used 
in NG is two times smaller than that of TIP3P model 
in earlier studies. Overall, for one helical turn, the effec- 
tive sampling volume in NG probably was several times 
smaller than before. For longer DNA the sampling was 
largely insufficient because the relaxation times for twist- 
ing and bending grow with the DNA length as L? and 
L 4 , respectively. For 5 turns this gives relaxation times 
comparable with and beyond 100 ns. 

The second note concerns comparisons with the WLR 
model. The authors analyzed MD trajectories with meth- 
ods that were never tested for WLR and there are strong 
reasons to believe that, in fact, the results would not be 
very different. Notably, the slowest stretching normal 
mode represents an accordion movement, with the whole 
WLR extending and shrinking concertedly. This mode 
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FIG. 1: Color online. Length dependence of the bend angle 
for Brownian dynamics of a WLR [jj. The same trajectory 
was analyzed with angle 9 measured between either true Z- 
vectors (closed red circles) or biased Z'-vectors (open blue 
circles) shown on the right. The bias angle <f> — 10°. The 
straight red line shows the theoretical dependence. 

can be extracted using the principal component analy- 
sis, and superposition of representative structures would 
give a pattern similar to that shown in NG, but it is 
not a sign of real correlations. The bending contribu- 
tion to the end-to-end distance is accumulated through 
the whole contour length and it is only partially elim- 
inated by using covariances measured between the ter- 
minal base pairs. Finally, oscillations with the helical 
period observed in some plots are inevitably produced 
by the algorithm used for constructing reference base- 
pair frames. This algorithm yields Cartesian frames that 
fit the global helical axis in an early ideal B-DNA model. 
For MD structures such frames are systematically biased. 
When this feature is taken into account for WLR it pro- 
duces periodical patterns similar to those shown in NG 
(see Fig. [lj. This is the main source of spurious oscil- 
lations in MD data . but it is an artifact rather than 
an evidence of correlations and deviations from the WLR 
model. 

In summary, all-atom MD simulations with the current 
AMBER forcefield somewhat overestimate the bending 
and twisting rigidities of the double helix. MD simula- 
tions agree well with the WLR model, with no detectable 
correlations beyond a few base pair steps. However, for 
non-integral number of helical turns the data analysis 
requires special care [!, 0]. Therefore, additional work 
is necessary for improving both the MD forcefields and 
interpretations of experimental data. 
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